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ABSTRACT 



The inherent symmetry of a spherical reflector microwave antenna 
suggests that if the fields incident upon the reflector from a feed 
required to produce a specified reflected field can be determined > 
then scanning of this reflected field nay be accomplished by rotation 
of the feed assembly rather than by rotation of the entire antenna, 

A technique has been determined to solve for the fields in vicinity 
of a hypothetical transmitting feed that will, follov/ing reflection 
from the surface of the spherical reflector, yield a specified field 
pattern at the reflector aperture. From given fields over a sphere, 
a portion of v/hich coincides v;ith the surface of the spherical reflec- 
tor, a spherical harmonic series expansion is used to calculate the 
fields on a smaller concentric sphere encompassing the feed region 
in solution of the boundary value problem. 

This technique is extended to a case where it is shown that 
geometrical optics approximations are valid over a specific region 
of the reflector system, and that by ray tracing, the fields on an 
arbitrary reference splicre, v/hose radius is confined only to the 
axis of symmetry of the reflector, can be determined * From this 
field distribution, the spherical harmonic series expansion can be 
used to detemine the fields on the surface of a smaller concentric 
sphere enclosing the hypothetical feed. It is shovm that this tech- 
nique has flexibility v/herein a number of origins could be examined 
for various feed spheres with fev;er terms of the series expansion 
than v;ith the fixed-origin case. 
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GENERAL BACKGROUND 

The fact that a reflecting paraboloid of revolution, when il- 
luminated by a source of spherical \?aves placed at the gconctrical 
focus of the paraboloid, vrill produce a uniform plane v7avc over the 
aperture of the reflector is vzoll kncv7iu The generation of a pencil- 
shaped beam antenna pattern has reached a state of refinement such 
that extremely narrov; bcamwidth, low sidelobe levels, reasonably 
broad-band operation, high aperture efficiency and gain are readily 
attainable for many applications. 

The axis of this pencil beam coincides with the axis of revolution 
of the reflector, and attempts to direct the antenna pattern by off- 
axis movement of the point-source feed have failed for all but the 
smallest off-axis angles. Therefore, to direct the pattern requires 
that the entire reflector and feed assembly be rotated or moved to aim 
the antenna pattern in the direction of interest. 

Widespread use of large aperture paraboloidal reflectors with wide- 
angle scanning requirements has stimulated considerable interest in 
development of improved feed systems to alleviate problems of ground 
noise radiation, spurious polarization effects and low aperture ef- 
ficiencies encountered in radio-tclcscopc applications. Other such 
problems that are encountered include difficulties in maintaining 
mechanical and structural tolerances of the reflecor surface at partic- 
ular frequencies of interest. 



The inherent syirjnetry of the spherical reflector suggests that it 
could provide some relief from the problems previously described. Assum- 
ing an energy source could be obtained that vrould, upon reflection from 
the spherical surface ^ provide an antenna pattern whose beam coincided 
V7lth the radius of the reflector upon which illumination occurred, it 
V70uld be no longer necessary to overcome the inertia of the entire 
system to meet a V7ide-angle scanning requirement, but merely to direct 
the feed assembly as desired. In such a case, it is feasible that large 
apertures could be illuminated to provide desired gain . Further, the 
symmetry of a spherical section enhances both ease of fabrication and 
verification of system tolerances in many applications. Other advantages 
that may be realised from such a system follow directly from speculation 
of possible application, includi.ng for example an antenna system whose 
reflector cou3-d be imbedded in the earth, V7lth few moving parts, pre- 
senting a low silhouette jto V7ind forces. 

Unfortunatly, analysis of the spherical reflector reveals its 

inherent spherical aberration. For a plane wave incident upon a 

spherical reflector, rays close to the axis of symmetry intersect at 

a point defined as the paraxial focus, exactly half-way between the 

center of curvature of the spherical section and the vertex. Other rays 

reflected from the surface intersect the axis at points closer to the 

vertex. This phenomenon may be observed in figure 1, where rays are 

constrained by Snell’s lav; to reflect V7Xth an angle equal to that of 

incidence, as measured from the surface normal at the point of interest. 

A more complete explanation of this reflection phenomenon will be givei\ 
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This suggests that rather than a distinct point focus there exists a 
fuy.zy focrJ region. Illumination of the reflector with a point source in 
the focal region reveals a diffraction pattern with wider beamwidth, 
higher sidclobes and much less gain than a paraboloidal reflector of 
identical aperture. It is this aberration that has presented an obstacle 
to development of the spherical rather than paraboloidal reflector for 
use as a micrcv7avc antennal Since introduction of the problem to the 
nicrov/avc community in WorliWar II, hox/cver, there has remained an 
acti.ve interest in the problem. Most recently, with construction in 
1960 of the Arccibo (Puerto Rico) Ionospheric Observatory’s 1000--foot 
diameter spherical reflector^ considerable interest has been generated 
in tVic solution to problems of spherical aberration at radio frequencies 
by means of various feed systems. 

1.2 It is appropriate to review past efforts relevant to the problem. Early 

attempts to explain optical pheiioncnn followed from geometrical theory. 

Euclid presented an analysis cf the reflection laws around 300 B.C, while 

observing that light energy could be concentrated by means cf 

S 
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reflectors* Early astronomers used spherical mirrors as 
collectors of light incident upon their telescopes « Presumably the 
ease of manufacture of these spherical mirrors led to their v'idespread 
use 5 even though the geometrical analysis of the reflecting properties 
of paraboloidal surfaces had been conducted* The difficulty of manufacture 
of the latter apparently precluded exploiting their point-focus charac- 
teristic# 

The concept of a spherical reflector for microwave use evolved from 
the World War II advances of radar technology* In 1941, Ashmead and 
Pippard suggested two approaches to the aberration encountered: 

1# For reflectors of large radius and restricted 
aperture, the surface of the spherical reflec- 
tor departs only slightly from that of a para- 
boloidal reflector# It follows that the diffrac- 
tion pattern is only slightly different from that 
of a paraboloidal reflector# If the sliglit depar- 
ture from plane x-jave illumination at the aperture 
and the value of gain accompanying the restric- 
tive aperture can be tolerated, this approach 
permits a simple design and is suitable for V7ide 
angles of scan# 
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2, For the general case of spherical reflectors p 
should the primary energy source be a special 
radiator that provides not spherical waves but 
waves v’hose departure from spherical waves is 
just that necessary to compensate for the differ- 
ence betv/een the sphere and paraboloid, then 
reflected wave fronts will be plane waves, uniform- 
ly illuminating the aperture* 

Early v/ork indicated that it was possible to create non-spherical 

v;avcs of approximate],y the right shape by proper combination of V7«ive- 

guide horns and dielectric lenses. Unf ortunatly , difficulties in 

Correcting for the larger phase error did not, at the time, justify 

pursuance of this approach. Rather, to meet the expedience of the 

v;artimc effort, the first approach of toleration of error rather than 

compensation was used exclusively. Subsequent studies exploiting this 

4 

technique r^cached a natural conclusion vjith the v;ork of Li , v;ho 
extended and verified the theory of restrictj.ve aperture large radius 
spherical reflectors. His experimental v7ork indicated that v/ith a ten- 
foot aperture hemisphere absolute gain comparable to that of a forty- 
inch paraboloid was obtained over a scan of 140^ v/ith rotation of 
the feed at a constant radius from the reflector. 

The second, more interesting approach to the solution allov7s the 
use of a larger aperture at the expense of increased complexity in 
cost, design and some reduction in the available angle of scan. 

io 



In tbo context of the letter approach corrective devices have fallen 
into the follcv/ing categories: 

1. lens structures 

2. auxilary reflectors 

3# line-source feeds 

4. transverse aperture feeds 

Lens structures j V7hile providing satisfactory correction eat optical 
frequencies j provide narginal correction in raicrovzavc applications. Sizes 
comparable to the reflector and primary source, requirements for large 
pieces of carefully controlled precision-cut dielectric and inherent 
dielectric losses, to say nothing of the costs and v/eights involved, 
relegate this category of corrective devices at microv/avc frequencies to 
a position of academic interest. 

Auxilary reflectors arc of minor interest \/heu the undesirable charac- 
teristics of shadowing or aperture blockage arc iimninent. Even in cases 
where this is of small concern, the excessive \7oight of the secondary 
reflector must be considered, as must the mechanical arrangements for 
support of the assembly, A gregorian reflector has been described by 
Holt and Bouchc^ v/hercin a reflecting surface intercepting a great deal 
of the energy incident from the spherical reflector is installed. The 
surface of the sub-reflector is such that it provides a constant path 
length from sene reference plane near the aperture to the primary source. 
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Experirasntal results indicate that the scanning capability v;as realised 
to 30^ by rotation of the sub-reflector on a constant radius froi^i the 
spherical surface. Sidelobe levels were very high in co:np>arison to 
those available with paraboloidal reflectors. V7hilc aperture blockage 
v;as estimated to be 1,5%, no figures on gain vjere given. 

The earliest atterapts at corapensation involved the fact that rays 
incident upon the spherical reflector were reflected to cross the axis 
of symmetry. It V7as postulated by Spencer, Sletteu and Ivalsh^ that the 
incoming energy could be recombined with an appropriate phase shift 
to provide a maximum of received signal. Reciprocity v;as invoked to 
suggest that a properly phased line source could produce a vjavefront 
that v’ould provide uniform illumination over the aperture. Techniques 
of geometrical optics V7ere used for path length calculations. The ~ 
sources considered were open-ended waveguides, poly rod travelling v7ave 
structures, discrete dipole sources alog the axis and slotted vzaveguide 
feeds. Indeed; the first feed system constructed at the Arecibo 
reflector was based upon a modification of this design! Aberrationless 
scan of over 60° was predicted for this initial installation-as V7cll 
as reasonable aperture efficiency. Actual rcsvilts were \;ell under those 
anticipated. The high-pov/er dual circularly polarized feed illiminated 

o 

the full aperture of 1000 feet with an efficiency of 21% at A30 MHz. 

9 

Love described a composite radiating system vzhcreln channel guide 
feed and a polyrod array provided a travelling v7avc structure to illumin- 
ate from both the line focus and the paraxial focus. He predicted a scan 
of 110° for a hemispherical reflector, and efficiency on the order of 



55%* Sche.ll‘^^±nvestigated the fields along the axis by application of 
diffraction theory. Ha dcmcnGtratccl the deviation from the gcopietrical 
optics solution in the non-zero wavelength and agreement in the zero- 
wavelength case. Schell solution consisted of the use of a sufficient 
nmiber of discrete sources located along the axis^ illtiminatlng partic- 
ular regions of the reflector. The analysis x;as based upon the fields 
present along the focal axis V7hen the reflector was illuminated by an 
incoming plane wave. In an attempt to solve the problem based upon 
Schell’s solution, Cohen and Perona^^ constructed a feed system at the 
Arecibo reflector* Their compound antenna feed consisted of a four- 
element linear array vrith a reflecting screen and a sixteen-element 
linear array located on the axis. Experimentation determined that gain 

V7as lovjer and sidelobes higher than had been anticipated. Aperture 

12 13 

efficiencies were on the order of 30-38% . McComick followed this 

study with an analysis that considered the transmitting case and, in 
particular, the aspect of polarization. His analysis concluded that a 
gain deterioration \jclo inevitable vlth a circular v;aveguide v/ith cir- 
cumferential slots that did not consider the azimuthal fields as t^ell 
as the longitudinal fields. The theory of McCormick v;as used by 
Love and Gustiiiic^^^ to illustrate the feasibility of illuminating from 
a leaky cylindrical waveguide that considered these fields. Host 

recent literature indicates that a line source feed v;ith high performance 

15 

has finally been achieved at Arecibo. LaLonde and Harris note that 
their feed has been installed at the observatory and has realized an 
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overall effixicncy of 70% at 338 MHzc Designing for a flat waveguide line 
feed v^ith transverse slots was necessary to avoid problems of previous 
feeds which j it v;as felt^ did not provide sufficient control over 
amplitude and phase variation along the line focus nor sufficient 
isolation from the effects of mutual coupling between the radiating 
elements of the feed* The lower frequency of operation v;as specified to 
avoid losses due to surface imperfections of the ixf lector. The design 
of the phase and amplitude required v;as based upon a geometrical optics 
analysis of the rays intersecting the axis. This design is significant 
in that it represents the first installation of a corrective feed in 
a spherical reflector that has yielded aperture efficiencies approaching 
those of paraboloidal reflectors . Aside from strengthening arguments for 
fixed-ref lector scanning antennas^ providing some 40 ^ of scan with high 
efficiency at Arecibo, the success of this feed indicates the potential 
of the spherical reflector antenna for fully steerable arrays j v;hero 
desired. 

It has been noted by Ricardi^^ that line source radiator correcting 
feed systems must be capable of producing a v/ave whose polarization and 
shape arc not easily obtained v;ith elementary radiators .This basic fact 
has caused present day line source correcting feeds to fall short of 
expectations, with exception of the LaLonde^-Harris feed. It should be 
noted, however, that the requirements of Arecibo call for high power 
with capability for dual circular polarizatiJr?. The line feed of LaLonde 
and Harris achieved success in providing for linear polarization. While 
the high power requirement may be easily met v/ith the slotted guide, 



the authors note that polarization flexibility has yet to be developed. 
Prccumably, the transverse aperture corrective feed designs have i\o 
inherent difficulty in polarization flexibility. Early nttcmptc to 
provide aperture corrective feeds mat with reasonable success in over- 
coming the aberration, but v7ere hampered by aperture blockage of the 
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inoderately large feed systems. Slctten and Mavroides v;orc succcsful 

in lovjering sidclobes substantially ^ v.^hile plagued ulth aperture 

blockage. An aperture feed v;as suggested for the spherical reflector by 

19 

Burrovzs and Ricardi in vzhicli they predicted a possible gain of 61 db. 

compared to the maximum gain of such an aperture the size of Arccibo at 

62*8 db. Their feed system vzas derived from a spherical harmonic expan-- 

sion of the fields necessary to produce a uniform distribution over the 

aperture j based upon the fields incident upon a spherical reflector. 

20 

Spencer and Hyde conducted a detailed study of the focal region fields 

associated with the spherical reflector illuminated by incoming plane 

vzaves and concluded that only a vector field solution considering polar-* 

ization vzould yield adequate information. They solved the diffraction 

integral by means of a stationar*y phase approximation , noting the 

similarities to the geometrical optics approximations. Their theory 

vzas shovm to have reasonable agreement with fields actually found in the 

focal region by experimentation, and their vzork has been used as a 

21 

basis for further study into aperture feeds . Their work, hovzever, did 

not offer a solution to the problem, but rather provided a theoretical 

analysis of the fields that could be expected. Work by Minnet and 
22 

MacAc Thomas has been extended to shovz that the fields in the image 
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of symmetrical focussing reflectors can be generated by certain hybrid 

23 

modes to provide theoretically aberrationless behavior ' . An aperture 
source is described that uses a corrugated vraveguide structure to generate 
a set of v;avcs that combines to provide a thcoreticall^y high efficiency* 
Further developing his earlier ideas of spherical harmonic expansions, 
Ricardi has recently provided a technique for synthesis of the fields 
on a surface enclosing a transmitting This treatraent for the 

transmitting feed synthesis provides for the solution to an electro- 
magnetic boundary value problem, identifying as boundaries the reflector 
surface and the surface of a smaller sphere enclosing the feed region* 

It is with this technique, substantially different from those involving 
geometrical optics or approximations of stationary phase, offering a 
solution that appears to contain inherent accuracy, v/ith \;hich this 
thesis is concerned* 
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CHAPTER II 



PROBLEM STATEl'lENT AND THE THEORY OF GEOMETRICAL OPTICS 

2.1 The solution proposed by Ricardi to the problcn of the spherical 
reflector nicrowave antenna relates tv7o sets of electromagnetic fields 
over portions of surfaces of t\;o concentric spheres. Due to the nathe- 
natical fornulntion of the problem^ one surface coincides ^?ith that of 
the reflector vrliile the other surface is taken at a radius of interest 
v/ithln that of the reflector. Origin of both spheres is at the center 
of curvature of the reflector (figure 2 ), The fields over the reflector 
surface are specified to be those fields that will produce a desired 
field distribution ever the aperture of the reflector follov7ing reflection 
from its surface. The fields at the reflector surface arc then expanded 
in a spherical wave expansion. Thus it nay be deternined that if tlic field 
distribution given by this expansion , evaluated at the inner surface ^ can 
be somehow generated by a source within this surface that the appropriate 
reflected fields v/ill appear, as desired, over the aperture. Interpreta-- 
tion of the fields at various radii has revealed that there are certain 
radii coincident v/ith nearly constant phase fronts. Synthesis of those 
fields V7ould be enhanced due to a potential lessening of the complexity 
of the radiating structure. Exactly which radius to be chosen for the 
inner surface is beyond the scope of this work, but is covered by 
Ricardi^^ 
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It is the hypothesis of this v;ork that the technique of solution 
to the boundary-value problem may be applied to a reference sphere 
vjhose orifp'.n is a variable distance along the axis of s^nnmetry and to 
the surface of v7h'Jch outgoing rays from a hypothetical feed located 
within this reference sphere would be traced if the spherical reflector 
were absents The r«ays outgoing from the region of the hypothetical feed 
would be precisely those rays representing the fields that v;ould> upon 
reaction v/ith the surface of the spherical reflector, reflect a desired 
illumination over the reflector aperture* These rays \7ould be based upon 
the rays reflecting from the spherical surface \A\en illuminated by an 
incident plane V7avc over the aperture. Strictly speaking, this analysis 
begins with the principles of the receiving case, but then requires 
the field distribution to be produced over the reference sphere by some 
transmitting feed. It is then the reaction of the radiated field with 
the spherical reflector that theoretically produces a desired v/avefront 
over the aperture* The geometry for this case is given by figure 3 . 

Thus, vzithin the liraits imposed on ray tracing technique by theory 
of geometrical optics, Ricardi*s results for the transmitting feed 
case may be compared to an approximate solution that originates v;ith 
the receiving situation. It is further possible that such a solution, 
knovm to be strictly valid only in the zcro-v7avGlength litait, may yield 
results v/ith less computation time than the exact solution within an 
appreciably small deviation in results. The choice of origin of the 
system of concentric spheres suggests that a more detailed analysis of 
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the fields in a particular region nay be conducted in an economical 
fashionj since the fields v?i3.1 be distributed over a sphere of smaller 
radius than the previous problem’s inner sphere v;ith origin at that 
of the reflector^ 

Uo thus proceed to describe the illumination of a spherical 
reflector by a uniform plane v;ave to illustrate certain character- 
is tics of the reflector « This is followed by development of geometrical 
optics principles essential to the description of the fields over 
the reference sphere. Finnlly^ expressions will be developed that 
yield the amplitude ^ phase , and direction of the electric field 
over the reference sphere. 



2.2 It has been noted previously that rays reflected from the spher- 
ical surface following illumination by a uniform plane v;avc over the 
aperture cross the axis of symmetry at varying distances from the 
paraxial focus. This phenomenon may be analytically formulated to 
illustrate the generation of caustic surfaces by these reflected rays. 
In fact, it V7ill be shown that the cavistic surface is the envelope 
of reflected rays. Consider the situation given by figure A. 
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Constrained only by the fact that tlie angles of reflection and incidence 
are identical , the equations of rays may be given by the slope intercept 
form of the straight line equation* That the rays are straight lines and 
the angler, equal v/ill be covered in detail in later sections* Viewing 
the system in the x-z plane v;e have: 

X « X 

tan 20 (2*1) 

z - z 

o 

Points on the reflector surface are given by: 

X ~ R sinO 
z = R cosO 

Letting points on the envelope be and substituting into 2.1: 
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( 2 . 2 ) 



X - R sinG - (R co 30 - z )tan20 
c c 

Thus points on the caustic surface as viev7ed in the x-z plane become 
functions of paraneter 0,. Wo then form the function f(x,z,0) ~ 0 that 
V7ill be satisfied by points on the caustic. 

fCxjZjO) = - R sin0 (R cos0 « z^)tan20 ~ 0 (2.3) 

Adjacent points on the caustic generated by rays incident at angle 
0 -{-60 will still satisfy the function for small 60. Since both functions 
f(x^, 2 ^^ 0 ) and f (x^,z^p0'l'60) are equal to zero^ their difference as v;ell 
is equal to zero. We can then form the expression valid in the limit of 
small 60: 



lim60*>O 



0-1-60) “ f(x^,z^,0) 
60 



= 0 



or 



d f(x^,z^,0) 



« 0 



(2.4) 



Differentiating and solving for z^ yields: 



~ R (2 sin^0 l)cos0 
2 



-2 2 . 



(2.5) 



This expression for z is then substituted into 2*2 to solve for x « 
Manipulation with trigonometric identities yields the parametric form 
of the equations for the caustic surface: 

X ~ R sin0 R( 3 sin0 - sin30) (2.6) 

c 

z “ ^cos ( 2 sin^0 1) = 3 cos0 - cos30) (3.7) 

^2 A 

Thus it has been shovm that the reflected rays’ envelope forms the 
caustic surface and quantitative expressions for this surface have been 
developed. They will be called upon in conjunction v;ith devclopnei'it of 
the expression for amplitude. 



2^3 Electromagnetic field problems may be formulated in tv7o contrasting 

fashions I follovjing from a choice of application of the Maxwell field 

equations. In one treatment, relationships may be derived that determine 

v;hat fields will arise from a prescribed set of sources. The field and 

vector Helmholtz equations are integrated by application of a vector 

Green’s theorem, and express the fields at an observation point as the 

sura of contributions from the sources distributed throughout a particular 

volume and from fields on the surface of the volume arising from sources 

26 

outside the volume. This rigorous treatment indicates that the field 
represents a flow of energy outward from the region of the sources, and 
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that these sources and fields satisfy certain assuinpt j.ons of continuity 

over the surface of interest, as required by application of the Green’s 

theorem* In the second treatment, the integral relations derived by the 

first method may be applied to yield expressions for the field vectors 

at a specified point, given the values of tlie fields over a surrounding 

equiphase surface, or wavefront, v/ithout direct reference to the sources 

27 

generating the v/avefrontso These integral expressions are the analytic 
formulation of the Huygens-^Fresnel principle, vjhich states that each po5_nt 
on the equiphase surface can be regarded as a secondary source of elec-- 
trie and magnetic current and charge v;hich in turn gives rise to a spher- 
ical v;avelet* The fields at a point of interest may be obtained by a 
superposition of these x^^avelcts, considering the phase difference upon 
arrival at the point of interest* Hence, energy flow is a v/ave phenomenon- 
one v;hich may be studied v;ithout reference to the sources of energy. As 
before, the rigor of the solution depends upon adherence of the given 
field distribution to several constraints as imposed by application of the 
Green’s theorem. 

V/hile the Huygens-Fresnel relation provides expressions for solutions 
to the v;ave equation that satisfy Maxvzell’s equations, complexities of 
integration of these expressions often render solutions impractical, or 
require approximations to be made that eventually lessen the rigor of the 
solution onl> after substantial manipulation. In these cases it becomes 
convenient to approach the propagation of vjaves frora the standpoint of 
geometrical optics, wherein a rigorous solution is sacrificed for a 
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simpler formulation that offers a more iiiumediatc solution, subject to 
several additional constraints on the problem. In geometrical optics, 
successive positions of ccjuipliase v/avefronts of the fields and an 
associated system of rays are related c 



In an arbitrary medium the wave field is characterized by both a ray 
velocity and v;avc velocity at every point, where the ray velocity is the 
velocity of energy propagation, directed as the ray passing through the 
point of interest, and the v;ave velocity is the rate of displacement of 
the wavefront in the direction normal to the surface* The energy propa- 
gates at the speed of light in vacuum, c. The wavefront propagates at the 
velocity V, a function of the medium* 

We proceed to describe the relations among wavefront and phase, con- 
sidering v7avefronts generated by some source propagating in a steady- 
state radiation* Let the v7avefront at time t^ be the surface L(x,y 
a surface of coUvStant phase == (u)/c)L^ relative to some point* The wave- 
front at some short time 6t latey is the surface L(x,y,z) = 6L* 

The phase difference between successive wavefronts is then given by 
(co/c)6L. Since the vjavc proceeds from one surface to the next in time 6t 
while the phase at any fixed position changes at the rate w, this phase 

difference must be equal to to6t. Further, if 6 is the distance betv7cen 

* sn 

the adjacent surfaces and if v is the wave velocity* then v5t == 6 . 

^ * sn 
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We then relate these expressions: 



(w/c)6L u)6t ~(oj6 )/v 

sn 

or 6L/6 = c/v f index of refraction of the mediuDi (9*8) 

sn 

Note thvat expression 2 # 8 is the directional derivative of a scalar L 

in the direction of the normal to the vjavefront. Ry definition of the 

28 

directional derivative of a scalar function 



d(J>/ds ^ T*V(J> 

where T is the unit vector indicating direction in which the derivative 
is taken* By definition of the vector dot product, v;ith 0 the angle 
between T and 

dij)/ds = = |t||V(}.| COS0 = |V 4 i|coa 0 

For our case the direction of interest in normal to the front, and 0 = 0 . 
Hence ; 



6L/6 = IvlI - c/v- n 

sn ‘ ‘ 



16? 



(2.9) 



Considering the curvature of rays in an iiihoiaogencous lUGcliurn we can 
demonstrate that the rays are rectilinear in a hemogeneous Tnedium* Let 
S be a unit vector in the direction of a ray at the paint of Int crest « 
This vector is norr:ial to the wavefront and h«as tlic direction of VL^ Then, 

|vl| = p or |vl| /n = 1 

Then S = VL /n (2.10) 

Let N be a unit vector in the direction of the radius of curvature of the 

ray at the same point of interest and p the radius of curvature of the 
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ray at the point. The vector curvature is defined' as N/p . This vector 
curvature is also given by dS/ds, where s is a distance measured along 
the ray. By definition and vector identity the following expansion can 
be made: 



dS/ds - (S*V)S = S X (V X S) - N/p (2.11) 

This cxpiression can be manipulated to chow a simple relationship bctv7c.cn 
curvature of the rays and index of refraction of the incdium. Dctail5; arc 
left to appendix 1. The result of intcre.st is: 

1/p = N*V(ln n) (2.12) 
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Thus for a homogeneous medium where the index of refraction is constant ^ 
or independent of position j v/e see that the right hand side of 1.12 
is zerop indicating that the radius of curvature of the ray is infinite, 
indicating that the rays are, in fact, rectilinear. It follovjs that 
VXS Oj \jhich is a sufficient condition for the existence of a family 
of surfaces orthogonal to a field of vectors S, 

The relationship of the amplitude of fields on successive VTavefronts 
can be deterrained from the premise of geometrical optics that rays are 
lines of flow energy and that no power flovjs across the sides of a tube 
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of rays. These results V7cre verified in the preceeding paragraph. Silver 
considers two surfaces of vjaves Lj^ and and a tube of rays that cuts 
out elements of area dAj^ and dA^ on the respective surfaces. By con- 
servation of energy requirements, the flow across any section normal to 

the tube will bo constant. In terms of the Poynting vector defined by 

% 

s = i^(EXir) 

V7C can formulate the flov; constraint by 
Si dAi “ dA^ 

Since in free space the Poynting vector is proportional to the magnitude 
squared of the field, between successive v/avefronts the follov/ing relation 
will hold: 



16 






(2-13) 



lEj^dAj = 



Considering the ray through point A on the curfiice coincident vritb 
the z^-axis in figure 5p let the planes coincide with the principle 

planes of at A« A ray through an adjacent point B lying in and the 
x-z plana \;ill intersect the ray through A at the point 0^, a distance 
ax'7ay« This distance thus is one of the principle radii of Lj^ at A, 
A similar analysis reveals the definition of the other radii R 2 based 
upon intersection of rays passing through adjacent points in the x-y 
plane. The point A’ on L 2 lies on the ray through point A, At A’ local 
coordinates x’, y’ are constructed. Since the fronts are spaced apart 
by the distance P, the two principle radii of curvature of the v;avefront 

L 2 at the point A* are P and R 2 P# 

% 

Considering the elements of area and dA^ around the points A, A’ 
bound by the curve C, C* on surfaces and respectively, v;e can 
write: 



dAi « /(xdy-ydx) 
dA = /,(x’ dy ’-y’dx' ) 

Z C- 



(2.14) 



From inspection of figure 5, and by the law of sines v?c can relate 
various distances: 



2.9 




By relating dx, dx' , dy, dy' and substituting into 1.14 \7c obtain a 
relationship between the area elements; 



dA„ = dA (R, . P) 



^' 1^2 



Substitut'-ion of into 2*13 ve obtain: 



(2.16) 




(P ^ + P) (^2 






For the case of a wavefront diverging from a region of sources this 
expression relates the amplitude at the wavefront farther from the 
source: 



3o 



(2.17) 



E, 



2 



E 



1 "/ 



(R, -i- H- pT 






1 



2 



Fropi the preceedlng development it v/ould app(?,ar that geometrical 
optics is based on the local behavior of the wavefront as a plane 
wave. Indeed, this is the basic tenet of the theory. The plane v;ave 
solution to the v;ave equation may be used to show that in the limit of 
decreasing vjavelength the geometrical optics solution satisfies the 
Maxwell fields equations. The vector solutions appropriate to describe 
linearly polarized electromagnetic fields that behave locally like 
plane waves have the foimi; 



The amplitude vectors may be complex, but phases must be independent of 
position to fit the plane wave solutions. Substituting these forms of 
solution into the homogeneous forms of Maxwell’s equations given by 



E « A(x,y,z) c 



j (ut“k^L(x,y,z)) 



(2.18) 



H -• B(x,y,z) e 



j (iot“k L(x,y,z)) 
o 



V X E -i- jcoviH 0 



(2.19) 



-jcJcE -!-V X H « 0 



SI 




f 



and expanding the resulting expressions ^ these may be solved for the 
vector amplitudes A and B : 



A ( VL X B) V X B) 

toe jti)C 

( 2 . 20 ) 

B VL X A) - _3^( V X A) 

WU jw;i 



where k is the free space propagation constant, k == w/pc/ , 
o 1 r I o *000 

Equations 2*20 may be solved to yield an expression in terms of A alone. 
Using the definition of k^ rmd the index of refraction 

n = c/v ^ (ey/e \i) ^ 

00 

we obtain the following expression: 



A - 



( VL X VL X A) _1 ( VL X V X A + V X VL X A) + 




V X V X A) 



( 2 . 21 ) 



Expanding 2.21 vith the vector identity 



VL X VL X A ^ L(VL‘A) - A(VL‘VL) 



VC obtain for A ; 



A « -1 { VL(A*VL) - AlVLp} + 1 ( VL x V x A V x VL x A) 



+ 1 V X V X A 

0 



( 2 . 22 ) 



A similar expression results for B, replacing A v/ith B in 2.22. If 
L and the first and second derivatives of A and B are finite, the last 
two terras of 2.22 are of the order 1/h^ and 1/h^ compared to the first 
terra. As the wavelength approaches zero, approaches infinity, and 
we are left V7ith the following expressions; 

A ;^1_ {VL(A«VL) - aIvlP' } 

(2.2.3) 

B = “1 {VL(B‘VL) - BlVLp } 

For these expressions to yield identity of amplitude of the vectors 
A and B \je see that the follov7ing must be true for the zero wavelength 
case; 

A*VL 0 



B*VL = 0 

These tv70 conditions note that A and B are transverse to the gradient of 
L, or that they both lie in a plane transverse to the direction of propa- 
gation. Furthermore, ciiiec 2.20 expresses B in terms of A V7S have; 
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e 



(2.24) 



B = /c/p{(VL 



X A) ~ 1 V X A } 



n 




V7e see that in the zero \7avelcnsth case the second tera vauisheo 




perpendicular to A as v;ell as transverse to the direction of propagationj 
since VL/n is a unit vextor in the direction of propagation. 

Thus the field vectors of geometrical optics posceSvS the properties 
of plane waves in v/hat inight be called the far~zone fields, indicating 
great distances from the source v/ith respect to wavelength. The condition 
for neglecting terras of the order 1/k^ and for short v/avclengths 

require that all associated derivatives are finite* In the neighborhood 
of geometrical focal points or caustic surfaces, the function L varies 
rapidly, and geometrical optics predicts infinite amplitudes , which 
must be rejected on non-physical grounds. At the region of geometrical 
shadov; the approximations begin to predict erroneous results due to 
the rapid variation of amplitude over distances small compared to 
wavelength. In these regions the simple methods do not vjork, and dif- 
fraction and scattering theory must be used. 

Further work has determined that the Huygcns-Frcsnel relation in 
the zero-wavelength case approaches the point-to-point amplitude relation- 
ship predicted by the geometrical optics concept of energy flow in tubes 
of rays. Agreement holds in the far-zone approximation, at distances 
greater than wavelength, and reveals the constraint on amplitude 



3 ^ 




# 




A Vs normal derivative that, besides being finite everywhere, 



1 

A an X 

This is satisfied in the limit of zero wavelength providing the frac- 
tional change in amplitude over distance equal to that of a V7avclmigth 
is small compared with unity. 



Since the basic tenet of geometrical optics requires that the V7ave-' 

front behave locally like an infinite plane wave, scattering phenomenon 

may bo interpreted geometrically, and appropriate lavzs of reflection 

derived from the field equations and boundary conditions* That these 
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results are valid can be developed from other principles of optics 
including Fermat’s principle and the law of the optical path, demonstrat- 
ing the properties of optical rays and the reflection phenomenon, or 
Snell’s laws* 

We consider the two-dimensional analysis of an infinite plane 
vzave iiicident upon an infinitely conducting ( that is to say perfectly 
reflecting) sheet* The geometry of the problem is described by figure 6* 
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"fl^UV'C 
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The incident fields v;ill have the forra; 



H - i h. ^°^®i 

i y 1 



( 2 . 25 ) 



E, = ^(i h.cosO. . ^ 
i ^ xi i z 1 X 



, . u » ^ \ “jk z cos0_, jk X sinO. 
•J- 1 h.siu0,) c o i • o i 






v;licre ^ is the v;avG impedance of free space, equal to (jj /c ) 

o o 

The reflected fields \7ill have the follox/inc^ forms 



„ _ . I *^jk z COS0 -f jk X sinO 

H - 1 h e o r o i 

r y r 



( 2 . 26 ) 



= ^(~i h COS0 i h sinG )e*^^’"o^ 
r^xr r zr r 



3(o 



Since nt 1 = 0 for all x ,,„<j y rl.o boundary condition on the tan,o„tlal 
=on.p„ne.,t of electric field re^uirca that (E, -IE) 

' i i'n ■“ U « 

Then, equating reflected and incident- f-tni i ^ 

na incident fields appropriately we obtain; 



0 . « 0 

1 r 



h, « h 
i r 



c . « e 
i r 



Thus we note that the following boundary 
face; 



conditions exist at the inter- 



(E. + E ) £= 0 

^ r tang ^ 

(E.) = fE 

i^ior ^\^ior 



C2.27) 



These 



expressions may be given in terms of 



a vector relationship; 



n X ( E^. + E^.) = G 



n ‘E « n .E 
^ r 



(2,28) 



cx=„i„o this vector relationship to detoiitine an expression for 
the iield reflected frot, a surface in ter.s of the ineident field and 
t>d turface netasal- as=u„i „3 that the field acts locally „s a plane ,,ave. 
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With tha vector form of the boundary conditions we form the cross product 
of the normal n and the first expression of 2.28c We then multi.ply the 
second of 2*28 v/ith the normal vector u and add the tvro expressions. 

This yields: 



n X E X n -r (n«E )n ” -n x E. x n 'h (n®E.)n 
r r 1 

Expanding by the vector identity 

A >: B X C « (A*C)B « (A*B)C 
we obtain the expression 

E - -E. + 2(n*E,)n (2.29) 

r 1 1 

This may be interpreted to express quantitatively the fact that the 
reflected field experiences reversal of its sense of tangential 
components in order to provide a zero tangential field at the surface 
of the ref lector. 



Thus we have demonstrated that the approximate methods of geometrical 
optics involving point”to--point transformations along rectilinear rays 
betv7ccn successive \-7avcf routs or surfaces of constant phase can provide 
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adequate phase aud amplitude information providiug certaiii constraints 



arc observed. Furthermore, that those fields satisfy Maxv^ell's fields 
equations and the wave equation within the limits of certain constraints 
has been demonstrated by a sero'^vravelongth limit. Ue su-mTnari^;c the results 
of geometrical optics in x;hich vc are j,ntcroctcd as relevant: to this \;ork5 



1. Energy passing through consecutive cross- 
sections of a tube of neighboring rays is 
constant* 



2« Rays in a homogeneous medium arc rectilinear. 

3. At a surface of reflection, an incident 
v;avefront behaves locally like a plane 
wave, and^ the reflection of rays may be 
geometrically interpreted in accordance 
v;ith specular reflection laws. 

The following constraints must be observed for the geometrical optics 
approximations to holdi* 

1. All lengths of interest, radii of curvature, 
for example, must be large compared to a 
v/avelength* 
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2« Spacing beUv^cen neighboring rays mur.t 
change very little over distances on the 
order of wavelength • 

3« The fractional change in anplituac of 
field quantities over distances on the 
order of wavelength must be small com*- 
pared to unity* 

4. Geometrical optics cannot be used in 
the region of geometrical cJiado'v/* 

5* Geometrical optics fails at focal points 
or caustic surfaces, predicting non-- 



physical results. 



2 ^ 4 It remains to utilize the previously developed expressions to describe 



the phescj, amplitude and direction of the electric field over the ref-- 
ercnce sphere « The choice of radius of the reference sphere is discussed ^ 

is the technique of ray tracing to be used* ExprosGioriS for phase 
and amplitude on the reference sphere are developed ^ relating these 
quantities to those incident upon the reflector surface « Finally ^ the 
distribution over the reference sphere is related to a coordinate system 
v;ith coincident origin « 

For the geometry of figure 7, we consider dctermiiiation of the fields 
produced by some source ^ directed towards the ref lector ,, that will produce 
illupiinatiou over the aperture in the form of waves propagating in the 
- z direction* The fields propagating in this direction v:lll be of the 
form: 



T- IT jh z 
E, H cc o 



or in spherical coordinates: 



E* n « o 



jk R COS0 



Since v;e desire to specify these fields over a reference sphere of radius 
R* v.’lth origin at z = cR, where c is some constant fixing the origin^ 
the phase expression must be corrected to account for the additional 



4 







rOUj 1 






I 







path lengtii of travel to the reference sphere. Inspection of figure 7 
indicates that the phase variation has tV 7 o cases of interest, one for 
angj.os 0 > 0 ^ other for the case of 0 < , v?hcrc ia the 

angle at v/hich the reference splicra intersects the spherical reflector. 
The variation in field expressions is introduced as? 



0 > 0c 



E^H « 



jk (Rcog 'I- L,) 
o 1 



0 < 0 ^ 



E,H « 



jk (Rcos 
o 




( 2 . 30 ) 



Inspection of the exponent suggests that if appropriate values of other 
pararaciters could be found so as to l:ccp the difference to a ininir.mm that 
synthesis might be enhanced, since the smaller exponent represents a 
slov 7 er varying sinusoid. Furthermore, amplitude conside'rations in the 
geometrical optics approximation, naraely the desire to avoid the region 
of the caustic surface, suggest that a V7ise choice v/ould be to consider 
the reference sphere intersecting the reflector surface at the latter’s 
extreme edge. Thus the radius of the reference sphere will be constrained 
by the radius of the reflector under investigeition, R, the position of 
the origin of the reference sphere z cR, and the maximura angle sub- 
tended by the reflector. By lax- 7 S of trigonometry x;e observe: 



r ’2 — 2cR^cos0 

m 



( 2 . 31 ) 
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The ray trac5.ng technique proceeds from a ray incident upon the 
reflector at a particular angle of interest. The ray reflects in ac-~ 
cordnnee v;ith the previously developed laus of reflection. VJe then 
assume that if a hypothetical source could produce the appropriate 
fields that vjould give the same illumination as vjc have upon the aper-* 
turoj these rays would remain as we have deterrained, but the sense of 
propagation of energy along the rays v7ould be reversed. Ue then remove 
the reflector, and trace the rays until they intersect the reference 
sphere. Within the limits of geometrical optics, vze can estimate the 
amplitude and phase at each point of intersection, based upon the 
extra path length travelled by the \7avcs, along the rays, to the refer- 
ence sphere. 

It is necessary to calculate the distance L. VJe proceed with an 
equation determining the distance betv;cen points on the reflector and 
points on the reference sphere, as measured along the ray of interest. 
Referring to figure 8: 

- (x - X (y - y -I (z -z (2.32) 

a m a m 

where 

X - 

a 

y ” 

z - 
a 



R’ sinO’ cos(J)’ 
R’ sinO’ sin;})* 



X = R sln0 cos6 
m ^ 



y - R slnO sind) 



R* cosO* -y cR 



m 



R cos9 



Noting from nycimetry that “ ({i', and substituting, 
wo obtain; 



expanding terms 



/' 



(2.33) 



(R'siuO* " P>.sinO)^ i- (R'eosO' •!- cR •* Rcos0 



To proceed, we must relate 0 and O' . Vievjing in the X"Z plane, the 
equation of ray is given by; 



X ~ (z - z ) tan20 
o 



y “ 0 



( 2 . 3/0 



From the equation of a spherical surface ; 



R'^ = -P (z - cR)^ -P y^ 



(2.35) 



Wc can solve 2.3A and 2.35 in the plane y = 0 to yield an expression 
represcutiiig the intersection of a ray and the surface of the reference 
sphere. Solving for intercept of the z-axis v;e obtain; 



R z 

sin (if- 20) sinO 



(2.36) 



Substituting 2.3A into 2.35, and expanding in a straightforward fashion 
we obtain, vfith substitution of 2.36; 

H5 



•• 2z( Rsin0s5.n26 -I- cRcos^26 ) •>' R^sin^O - cos^20(R'^ - c^R^) - 0 



DctaiXn of the rannipulntion ore included in appendix 1. Nomalizing a3.1 
distances to R, radius of the reflector imdcr investigations v;e. obta5.n: 



z 2z ( sin0siu20 *1 c cos^-20 ) -’r sin^O- cos^20(R'^ 
n n n 






c^) - 0 



(2.37) 



v?here 



R' = R' 



z = z 
n — 



Solving the quadratic 2,37 by usual methods, \;a obtain for z^ : 



z = ( sin0sin20 c cos^20) -1 
n 



(siiiOsin20 c cos^20) - sin^Of cos“2C(R^*^^ - c^) 



(2.38) 



Noting from the geometry of figure 8: 



R cos - Z 2 - cR 



COS0‘ = 



(z - cR) 



“ ~ 

R' 

n 



HG 



or 



(2.39) 



Substituting equation. 2* 38 into 2«39 the follovring rclaticnship batv/eeu 
angles 6 and G* is obtained: 

coc6* - ( sinCsin20 d- c cos^20 - c 

u 

+ ^(sin0sin28 *h c cos^20 ) - sin^ G -h coc^-20 - c^) ) 

(2.40) 

Equation 2.40 states that there is a relationship betvjeen angles such 
that for every angle of incidence on the reflector surface there exists 
an angle G’, as measured from the origin of the reference sphcrCj that 
describes the intersection of the reflected ray traced back to the refer- 
ence sphere and the surface of this sphere. Thus the angles may be 
related in terms of the parameters of the reflector under investigation^ 

R and 0^ , and the center of the reference sphere implied by the. value 
of c. Once the angles associated x^rith a particular incid ent ray have been 
related, they may be substituted into 2.33 to yield the distance L as 
a function of the point of interest on the reference sphere. 

Equation 2.17 \*7as developed to express the amplitude of field quanti- 
ties on one vjavefront v/ith respect to the field quantities on another 
front. It is seen that the amplitude varies V7lth the tv7o principle radii 
of the x;avefront and the distance of separation between fronts. It is 
necessary, then, to develop expressions for the radii. 




48 



Onca again v;a consider the ease of a spherical reflector il3.uminatcd 
by an inconing plane vjavec A typical reflected vravofront lias been 
constructed in figure 9j vjhere the front represents a surface of 
constant phase as defined by the total path length measured along 
each ray, -I- ^ wavefront is formed by connecting normals 

erected on each ray at the point v/here total path length equals the 
radius of the reflector. Figure 9 shov;s the front as viewed in the 
plane y 0* From symmetry , it may be seen that the surface of the 
front can be generated by revolution of the projection in figure 9 
around the axis. 

At any point of interest, a v/avefront may be passed through this 
point by choice of the appropriate path length « At each point the 
v^avefront will have tv;o principle radii of curvature. Since a ray ±s 
normal to the front and any tv?o adjacent rays intersect at the caustic 
surface, forming the envelope that is the caustic surface, it is seen 
that one radius of interest is the distance from the front at the point 
of interest to the point on the caustic, as measured along the ray. The 
second radius of curvature may be determined from the axial symmetry of 
the system. In the direction, adjacent rays v;ill intersect at the axis 
of symmetry, and the second radius of curvature will be the distance 
from the front at the point of interest to the axis, measured along the 
ray. 
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Rays incident upon the reflccto’J nt an angle o nre reflected at this 
angle, pass tangent to the caustic surface and cross the axis. The 
coordinates of the point to vdiich tangency occurs arc; given by ec|uations 
2,6 and 2.7, the parametric equations of the caustic. Frcn figure 9 vje 
see: 



R. = MU = 

1 me w c 



(2.41) 



R„ = MVI = /(z T“F’T~lx - 0)‘ 

2 no m 



(2.42) 



where 



z - 
m 



X == 



m 



RcosO 

Rsin0 



z = R 



sin0 



siii2o 



Substituting into equations 2.41 and 2.42, and normalizing all quantities 
to the radius R we obtain: 



R, - 

In 



cosO 



(2.43) 



R 



2n 



2 COS0 



(2.44) 



Details of the maaipulaticn involved are covered in appendix 1. 



So 



Although tliG \7avcfront has been analyzed by considering illumination of 

tlic reflector by an incoming wave, it is the fact that if a source can 

be provided that can generate such a v;avcfront as has been describedj the 

wavefront \7ill be reflected from the spherical surface to yield an 

outgoiiig plane v/ave ever the reflector aperturOo Figure 9 represents a 

graphical illustration of the method of compensation suggested by 
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Ashmead aind Pippard, a wavefront vjhose departure from -sphericity is 
just that needed to transform into a plane v;ave upon rcflcctione 

Since we v;ish to describe the fields over the reference sphere as 
a function of coordinates with origin at that of the reference sphere, 
wc must introduce a transfonnation of the field vectors^ Initially, 
field quantities v^erc described upon incidence in terms of coordinates 
with origin at that of the reflectoTe The geometry of interest is 
described in figures 10 and 11, considerdng both cases where the center 
of curvature of the reference sphere is exterior and interior to the 
intersection of the reflected marginal ray (edge of reflector being the 
point of incidence) and the axis. Note that due to axial symmetry there 
is no modification necessary to the (}> components Sy inspection: 



i^,f « i^ cos (0^ 0) + ig sin (0* - G) 

1 I ^ - i sin (0* - 0) i_ cos(0* - 0) 
u 1* u 



(2.45) 
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Solving for the unpriraeci components in terms of the prined cystcra: 



i^. ~ i^_,cos(0* - 0) “ i^,sin(0* “ 0) 
ig = i^.,sin(0' G) i^tcos(0’ - 0) 



(2.46) 



From the preceeding derivations it may be seen that equations 
2#17j 2c30 and 2.29 provide sufficient information to form an expression 
for the geometrical optics amplitude, phase and direction of the electric 
field distribution over the reference sphere, based upon a desired 
reflected field . Quantities of interest, including the radii of curv-- 
aturo of the v;avcfront, distance from the point of the reference sphere 
under consideration to the corresponding point on the reflector, and the 
relationship between angles are given by 2.43, 2.44, 2.33 and 2.40 
respectively. Final3y, transformation of the vector components given by 
2.46 relates this field distribvition to coordinates at the center of the 
reference sphere. Thus the geometrical optics field is described over the 
reference sphere as a function of the i*eflectcr surface under inveotiga^ 
tion and the particular origin of interest of the reference sphere. 
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In summary, we formulate an expression for the general case i/here 
a desired electric vector field is specified over the aperture c Choosing 
this field to have imity rmiplitude, V7e describe the field that v;ould 
occur over the refereuce spheres 



Amplitude: E 



reference 'reflector 
sphere 



R, R. 

In 2n 



'^( 2 . 17 ) 



(R, <1^ L ) (R^ L ) 

in u 2n n 



Phase: 



r»u / \ jR R( COS0 L ) 

Phase ( E - ) « o ^ 

reference 

sphere 



(2.30) 



Direction: E- « *-E*h2(n«E)n 

reference r r 

sphere 



(2.29) 



Combining: 



E . 

rercrence 

sphere 



Rs R^ 
In /n 



(R- h ) (R^ L ) 



In 



n 



2n 



1k R(cos 
o 



{~E -!• 2(n*E )n} 

r r 



(2.A7) 



S‘i 



v/hcre; 



COS0* - 




cos6 

2 



(2.43) 



R 



2n 



.J 

2cos0 



(2.44) 



(2.31) 



L = (sin0 - R’sinO’) •!- (R' coc0' + c *- cosO) (2.33) 
n n n 



1 ( sin0siii20 -!• c cos^20 

R' 

n 



c 



r\ 



siu0sin2O + c cos^26)^~ sin^0+ cos^-20(R^^ 



'?■) ) 



( 2 . 40 ) 



ij. “ i^icos(0' - 0) - ig,sin(0’ - 0) 
ig i^,sin<'0* " 0) ig,coG(0’ -• 0) 



(2.46) 
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2*5 In this section, v:e apply the results of the general expression 2«47 
to a particular case of interest* Specifying the field at the aperture 
shown in figure 12 to be a uniform plane wave linearly polarized in the 
X- direction V7c have: 



E 

aperture 



i 

X 



jk z 
c^ o 



Expanding in spherical coordinates: 



E 



aperture 



jk RcosO 
e'^ o 



i sin0cos<j) 4* i^cos9cos({> 
r 0 






(2. AS) 



By applying equation 2*29^ reversing the sense of tangential components, 
wo can determine the sense of direction of the field over the reference 



5^6 



sphere: 



dircction( E 



- ) - i siiiOcoGd) - l-cosOconA -’r i nind) (2.49) 

reference r ^ 6 ^ 

Siphere 



Applying equations 2.46 to 2.49 the direction of the field vector 

can be expressed in terms of the origin at that of the reference spheres 



ij. i^,cos(0' - 0) - i^,sin(0' - 0) 
ig “ i^,sin(0' - 0) + ig,cos(0' ^ 0) 



(2.A6) 



Substituting s 



direction( E ^ ) == i ,cos(J) ( sinOcos(0* G) - cos0sin(0’ - 0)) 

reference v v \ // 

sphere 



- i ,coG(J)( sin0sin(0* ^ 0) -y cos0cos(0’ - 0)) 
0 



•i- i siud) 



(2.50) 



Collecting terms and applying trigonometric identities for the difference 
of angles we obtain: 



direction 



( E 



reference 

sphere 



) 



= i^, cos(J)Sin(20 

-h i^sind) 



0 ') 



i ,cos(f»cos(20 - 6’) 
0 

(2.51) 



i 




Thus fron, previously Jovolo;.e,-1 rolr.ticr.ships 2.J7 and 2.30. oxprossions 
for .'.uplitudo, phase, and ditoctioi, of the oloctrlc field over the ref- 
ereaco sphere can be conblnej to yield the vector field solely as n 
function of the particular reflector under investi-ation and the choice 
of origin of the reference cplicre o-lcnf; the c^y.xs, 

Hlille sinplification could be pertorned on the realutinc oxpresoioa, 
the field over the reference sphere is. in itself, only a neans to 
eventually descrihe the fields over a smaller sphere surrounding sone 
source region, concentric nith the reference sphere. Since the sinpUfica- 

tlon would not enhance this procedure, it will „ot be pursued. Rather, 
the resulting expression 



reference 

sphere 



In 2n 



n 2n n 



y. 



■2 •• L ) 

O O II 



( i^,cos.fsln(2e - 9') - i ,cosdcos(29 - O') + 1 sInd, ) 

<P 

(2.52) 

(w.ieie relationships arc defined by 2,/)3^ 

2.44, 2.33, and 2.40 and are listed for convenience in section 2. A) 

itself readily to tecniques of nachine conputation. As a result, 
the following pages contain graphs representing the field component.- 
In anplitude and phase for a particular choice of paraneters. Th. 
r.re calculated for nornalized annlitude a.nd reflect no variaf^ 



S8 



ajr^inuLb.al 6 dircction« Thus the conpcnents E / and E. indicated ±n the 

U (J) 

gr^iph nur>t be rmltiplied by costj) and sinej) respectively to represent the 
true nature of the field over the entire sphere# Similarly ^ the E^, ccra- 
ponent must be multiplied by the variation cose}) # The computer program 
developed to calculate the field over the reference sphere is included 
in appendix II# 
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CH/iPTER III 



ELECTROMAGNETIC FIELD EQUATIONS, EXP/iNSION^ AND SOLUTION 

3*1 In the prececdxng chapter expressions for the electric field over the 
reference sphere x^erc deterrr.incdc In this chapter those fields x*;!!! be 
raatched at the reference sphere to an expansion of spherical hn.rmonics 
follovring from solution to the vector x^ave equation in spherical coordi-- 
nates. These harmonics x-jill be evaluated at a particular radius xvithin 
that of the reference sphere to yield the electric field over the surface 
of that smaller sphere enclosing some source region* These fields Xv^ill be 
compared to the solution of Ricardi for the same region* Finally, these 
results and the comparison of solutions Xv'ill be interpreted* 



3 A 

It may be shoxni that solution to the scalar Helmholt 22 equation 



( 1^2 )y 0 



(3.1) 



for spherical coordinates is given by separation of variables to be: 



m = b: (kr) P^'(cosG) 

ran n n sinmcj) 



(3*2) 



where bt^(kr) is the spherical Bessel function and P^^|(cos0) is the 
associated Legendre function of the first kind* The choice of cosm^i 
or sinm({> depends upon the variation of the fields at the boundary of 
interest. 



bZ 




Il 

Cf. 




Tha ficldr, under InvcstD gaticn in chapter 2 v;ere shown to be vector 
fields, and the solution to the vector HelnholtJi equation imust be 
considered as: 



(V^ k^)A ^ 0 



(3.3) 



where A is a vector fields representing either the magnetic vector 

potential or its dual^ the electric vector potentialc It has been 
35 

shown that a general solution may be formulated by considering a set 

of fields, botli magnetic and electric, vrhich are transverse to the 

radial vector, i . In this case the solutions to 3.3 are found to he: 
^ r 



A ~ r i 

r 



F = ¥ r i 



(3.4) 



vjhere is the solution to the scalar Helmholtz equation, with the superb- 
script indicating that is solution to the vector Helmholtz equation 
for the magnetic potential A or the electric potential F. 



G3 



The most general r;oluticn for fields in a source-free lioniogeneous region 

36 

in terms of vector potentials is found to be a linear combination of 
fields \-7ith the bcliavior of equations 3.4, v;here; 



E - -V X F - jojy A -f _ V (V • 

jwe 



H == V X A - jtocF + 1 V (V* F) 

Jup 



(3.5) 



Since the solutions for A and F involve the product of r and Y , it has 

been found convenient to adopt .an alternative definition of the spherical 
37 

Bessel function : 

H (kr) - kr h (kr) 
n n 

If) 

The superscript indic.atcs wc h.nve chosen tlie spherical Hankel function 

of the second kind, whose behavior' indicates prop.agation outvjard from the 
origin, .as .ipplic.able to this problem. Anticipating further, Ej.ncc the 
fields we wisli to expand are of the form of equation 2.52, v7horc the var- 
iation v;ith (j, is given by nini}) and coS(t> , V7e clicose n = 1 , Thus: 



F “ 
r 



E HJUr) 
n 



P (cos0) 



n 



{ sin^’ }. 
{ cos } 



n 



n n 



(3.6) 



L b. H (kr) 



l-hcoae)' ! 

n { coS(J } 




4 



Expanclinp, equation 34 5 in component form v;e obtain; 



“ „L. 

jUG 



A 

) r 







rsiii6 8(|i jwer 8r80 



r 


8F 

r 

80 


+ 1 

jwersinG 


s2a 

r 

8r8<^ 




82 

( ““ 


■y R2 f 

“■ ) r 





jtop 8 1-2 




r 8G jeoyrninO 8r8<J> 




rsinO 8({) jeoyr 8r86 



(3.7) 



Referring to the componoiit E„ in 3.7, we note that (j> variation for 
0F 

A and ^ ' r is the sane. Hence \7c choose the appropriate variation in 
^ 9 <{) ^ 

equation 3,6. Ue note further that the function P^(cos0) is identically 



zero for n - 0« The series need begin for n = 1, and nay be written: 



~ A (?'» , 

A “ E a II. (kr) P (cgg6) cos<{) 

r , K n n 

n=l 

CO ^ (2 ' 

F « E b II (kr) P^(cor>0) slnA 

r , H n n 

n-1 



(3«8) 



Substituting into the appropriate conpenents of 3«7 the electric field 
can be formulated : 



E 

n 



r. n(n*^'l) 
11 

jojcr^ 



A(a) 

H (kr) 
n 



11 



(cos0) 



COS({) 



i 

r 



b H (kr) P'(cor0) cos(J) i 
11 n u ^ i 



rsin0 



(3.9) 



MZ1 



jcuEr 



a cosd) n 
n — — 

dr 



1 



dll (kr) dP (cog9) 



n 



de 



^ JT>1 / rs\ 

V Tt /I \ t a <1^ (COG0) 

b H (kr) sinc^ n 
n n “ — 

do 



^ a siiKj) P^(cos0) 



jwer sinO 



dH (kr) 



n 



dr 






^il 










t 







For convenience \re define: 



A (21 




A (21 


d 11 (kr) 
nj 




k 


dr 




d(kr) 



k H (kr) 
n 



and since v;e are concerned v;ith the case of 



propagation in free space 



k ~ k =w/g y 
o o o 




Substituting the relationships from above, multiplying through by j/j - 1, 
and including the suTamation we form equation 3^10 from 3o9 : 



E 




n(n*M) 

n 






jtoer 



AC2’ , 

H (kr) P (cos 6) cos* i 
n n r 



( 3 . 10 ) 



+ H^(kr) 



, , F (cosG) .... r>J-/ . 

-cogO Cl n + siri(J) P (cog6) i 

Tq~— e n ^ 



A [Z) 

+ b H (kr) 



n n 






. T^l/ ON • t j A 1 ^ (cosO), 
•cos cj) P^^(cosG) sinc^ 

siu8 ^ d 0 ^ 



J 
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Introducing tangential spherical harmonics as 



defined by Morse and 



Fcshback 
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cosmej) n 

do 



dP (cos0) 




sinO 



Oai) 



c 



o 




Substituting 3oJl into 3*10 v;e can fern an expression for the tangential 

ficldvS. By exploiting certain orthogonality properties of the tangential 

harnonics. \7c can derive the constants a and b in 3c 10 from the fields 
^ n n 

specified at the bovmdary of interest. In our problem, it V7ill be the 

tangential fields on the reference sphere that v;ill be specified. Note 

that from equation 3.10 the two constants a and b are sufficient to 

n n 

describe the field vector. It must also be noted that development of all- 
expressions in this chapter has proceeded for the case where the origin 
of the sphere for vdiich boundary conditions vrere given coincided v;ith 
the origin of the coordinate systera for the spherical harmonics. It now 
bccomcG necessary to stipulate that the harmonics must be expanded about 
an origin at the center of the reference sphere. To adapt these e::pression 
\7e need only to substitute 0’ for 0, recogviising that the angle is do*“ 
fined as in chapter 2, and for the radius variable r v/e substitute r’, 
indicating that it is measured from the prime coordinate origin. 
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From the pro.ceeding vjc can write for the tangential field components: 



tang 



1 

r 



AC?-' 






T. /n(n-'rf) b H (kr) C® d- iCa H (kr) B, 

* n n In ' n n In 



n=l 



J 



(3.12) 



Coefficients £ira evaluated by forninn the vector dot product C_ *E , 

Iq tang’ 

multiplying each side of 3,12 by sinOdGdcJ), integrating over the 

interval 0-tt end 0-2tt on cj) rmd 0 respectively ^ evaluating the fields at 

the reference sphere^ radius P.* , The following properties are listed 
39 

below: 



/ / 



e e 

.B, 

Iq In 

c° .c” 

Iq 1ft 



sinGdGdcJ) 



2rn(n-!-l) n.-q 

2n-fl 



0 ftrq 



/ / ^in*^3q 



0 for all n^q 



Applying these properties to equation 3^12 we have; 
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i(2rL-!-l) R' // E sinO 'd0 'd<J> 

~ rt In 

C2Trn^(kR’) (n(n-'rl))^^^ 

(3.13) 

~(2u-!-l) R’ // E *C^ sinO'd0'd({) 

?,TrH^(kR') (n(u+l))^^^ 



where is the vector representation of the tangential ccvnponants cf 
the electric field over the reference sphere at radius R' e To perf orn 
the integration indicated in equations 3*13 V7c proceed by expanding 
the field in a fouricr series* We choose a series appropriate for 
the fields over the reference sphere, in this case a scries with even 
syr^imetry about O' - 0, and \7here the <p dependence matches that of the 
fields at the boundary, as suggested by equation 2*52* Proceeding: 

CO 

E (R’»0*>^) ~ r ( c cos(p i , d sind i.)cosni0' (3*14) 
rt ^ ' m ^6’ m <P 



The fouricr coefficients arc obtained in the normal fashion, and arc 
covered in appendix 2* 

Consistent V7ith the notation of Ricardi we define coupling coefficients 



/^cosmO' P^(cosG') dO* 
nn o n 



== /^cosmO' sinO' 
mu o 



d P^(cos0’) 
n 

do' 



dO' 



(3.15) 
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SubstitutioTi of 3.15 into 3.13, u?jing the definition of fouricr co- 
efficients from 3«14 va. obtain: 



a 

n 



H (kR*) 
n 



E 

n-O 



c 

Tfl 



nn 



d 3:^ 
la mn 



) 



b 

n 



- A R’ 

H (kR*) 
n 



E 

Ta=0 



C 1 
in mn 



d 

n 



mn 



) 



(3.16) 



where 



2n + 1 

2(n(n+l))2 



Therefore to fully express the electric field in the spherical harmonic 
expansion we need o'aly to evaluate the fouricr coefficients in 3.3 4 and 
the coupling coefficients given by 3.15. These coefficients, vjhcn sub- 
stituted into the scries 3.15 yield coefficients of the spherical 
harmonic expansion. 

These constants are normal3.y evaluated by a recursion tcchnique^*^ 
rather than solving for each term by evaluation of the integrals. This 
recursion technique is particularly easy to realize v/ith machine compute 
tioa, and is inherent in Ricardi*s method of evaluation. Thus coefficicn 
\/ill be computed rather than analytically determined. 

II 



In suTiimary the relationr.hips that express the electric field as a 
spherical expansion of a specified field on a given refercvicc sphere 
such as shoun in figure 12 are listed: 



K(r\e’,(i.) - 



(3.17) 



-I 


CO 




I. 


r' 


u-1 


r 


A ( 2 *^ 


b 


H (1 


L " 


n 



•a n(n-!-l) ,, 1, 

n H^^(kr') P^^^(cosO’) 

jwer' 



In 



jC a H (kr') B, 
n n In 



o e 

where and are given by equation 3*11 \7ith m-1 : 




cosjS^ (cos 0 ' ) 

sin0’ 




silKj) 



dP^(cos0 ’ ) 
n 

do' 



J 

/ n(lVrl) 



cosef) 



dP^(cos0 ' ) 

_ n 

cfo'" 



^ 0 » 




and and given by equations 3*16: 



a 

n 



j R* 

ji 



E 

n=0 



( c„ I 



,-d 



in rnn 



d J ) 

El tm 



7Z 




t 



00 



b = 
n 




II' 



riur') 



E 

ni-0 



c I 
ni ur 



d 

n 



nu 



) 



. ?.n 1 

== 

2(n(u+l))?- 



and the coupling coefficients 



mn 






„d . , 

\an 



3.15: 



“ /cocraO' P^(cos6') dO' 
rrm o n 



I 

nn 



/conn0 ' sinO * 
o 



dP^(cone') 

u 

de' 



de' 



The fourior expansion cosf f icent:s c and d are defined in equation 3.14 

rn m 



E , ^ E ^(R',0',4>) “ ^ ( c co£) 6 i., + d siiKji i,)co 2 inG' 

rt rt ' w ^0 n 

n-0 



n 



I 




E is the vector representation of the electric field tangential to the 
rt 

reference sphere ^ as determined from 2<52» 

While the series representation gives an accurate solution for the 

sum of all terns 5 it is recognised that this is impractical to realise. 

An important section of Ricardi’s thesis dealt with justification of an 

4] 

appropriate truncation of the series solution . That analysis deals v/ith 

subjects beyond the scope of the present work, although it v;ill be noted 

that the series truncation stems from a desire to avoid generation of 

highly reactive waves that do not, in practice, yield that much accuracy 

to the end result. On the basis of a rigorous analysis regaeding genera- 

42 

tion of these highly reactive v;aves , it will suffice here to state that 
satisfactory results are obtained with the spherical expansion v/hen the 
number of terns of the series is on the order of N ~ kR| , v;hcrc is 
the particular radius at v/hich evaluation is desired. 

Therefore, equation 3.17, evaluated at r’ = , v/here is the radius 

enclosing the transmitting feed v/hich will produce the desired field 



^reference 

sphere 

specified by equation 2.52, V7hich will in turn produce after reflection 
a desired field distribution at the aperture, yields 



feed - ■’>f. 6'. 

sphere 



The equations for magnetic fields over the feed splierc can be obtained 
froia Maxwell's equations* In particular: 

V X , 

feed — feed 

sphere jtop sphere 

In this analysis v/c have concentrated on solution of the electric fields 

at various boundaries. To produce the desired fields at the aperture it 

is not necessary to separately excite both $ since 

sphere sphere 

through Mnx\7ell*s equations ^ the existence of one field specifies the 
other. 



15 



3.2 Xo illustra.te th.e principles, outlined in preceeding 
sections, and to show that tlies'e. principles provide an 
economical and efficient method of mathematical expansion 
of fields from the geometrical optics reference sphere to 
an interior feed sphere, a particular case was chosen for 
evaluation , 

The parameters of interest were chosen for c = .52, 

representing a displacement of the origin of the feed and 
reference spheres to z = .52R for a reflector of radius 

E = 93X and subtended angle of 37,4 degrees. The field 
distribution of magnitude and phase were presented in 
graphical form in section 2.5 for this choice of parameters. 
Figure 3 of section 2.1 illustrates the geometry of the__^ 
eccentric coordinate system. 

It must be noted that the choice of eccentric origin 
chosen for this case of investigation was rather arbitrary, 
as was the choice of inner sphere radius. This choice was 
guided, however, from principles of geometrical optics. The 
rays incident upon the extreme edge of the reflector surface 
are reflected into the focal region, forming the caustic 
surface. The points at which the edge rays are tangent to 
the caustic surface define the marginal surface focus, that 
plane surface generated by the rotation of of figure 



15 about the axis of revolution of the reflector. 



Note that all reflected rays must pass through this surface ^ The 
smallest such surface through which all reflected rays must pass is 
defined by the circle of least confusiouj generated by rotation of 
A-A' of figure 15 about the axis of symmetry^ 

Wliile geometrical optics v;rongly predicts infinite fields in the 
vicinity of caustics, the fact remains that converging rays indicate 
an area of energy concentration. Hence, an intuitive choice of probable 
location for a minimum-sized feed v/ould be that aperture which most 
closely coincides with the circle of least confusion associated with 
a particular reflector. It is a simple matter to show from the geometry 
that a choice of origin betv^een z = «5R and *6R requires a feed sphere 
on the order of 5--10 X for a reflector of 93 X v;ith subtended angle 
37.4 as considered in this study, to meet these specifications. 
Hence, a representative distribution of fields over such a feed sphere 
that V70uld be required by the^ tangential boundary conditions to yield 
a uniform plane Vv^ave over the reflector aperture is given by figures 
16. Amplitude and phase of three components of electric field over 
the surface of a feed sphere of radius 8.4 X V7ith origin at z « .52R 
for a reflector of radius Pv ~ 93 X are presented in figures 16 A-F. 

This solution is a graphical result, obtained from evaluation of the 
series solution given by equation 3.16 and 3.17 of section 3.1. 
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3.3 For tho ch^oica of pai;ainaters indicated in section 3,2 
the magnitude and phase of the field distribution over the 
reference sphere, as given by equation 2.52, were presented 
in graphical form as figures 13 and 14 of section 2.5. The 
nearly identical distribution of tangential components and 
the near--vanishing radial component of the field are due 
to the relationship betvzeen angular variables 0 and 0*, 
where to first order ; 

20 = 0 ’ 

Only the tangential field components are used in the 
boundary condition to generate the fourier coefficients, 
as noted in appendix II. The tangential fields were expanded 
in the following manner -from equation 2.52: 



tan 



= ( 



Eecos^e 



E s ind) i ^ 
<)) (j) 



)e 






= E - CO s d) ( c o s ij; + 
tan 0 



jsin]|;) ig + 



E,sin(ti(cosii 



+ jsinif))i 



<!> 



S5 







A 



A 



P 



Tl\e. re.Al a,ncl injagin^j;y pa.rt:s, le.ss tiue cos(j)^ siu(p dependence 
were then plotted graphically and are presented as figures 
17 and 18 , Due to the similarity of magnitudes resulting 

from the near ly-ident ical nature of tangential components 
indicated by figures 13 and 14 of section 2.5, only two 
rather than four plots are presented here. The calculations, 
however, retain the existence of four separate components, 
two real and two imaginary, of tangential fields, These 
components were then expanded into a fourier series of the 
f orm : 



where, c , d are now complex fourier coefficients determined 



as noted in appendix II, and M is the limiting value of 
the highest harmonic desired, given, as noted previously. 



H 




m 



s in<t> i^ ) cosmG 



m m 



by 



M 



kr 



with r the radius of observation. 
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To check on the accuracy of the fourler expansion technique , the series 

% 

was suntmed using coefficients determined by the geometrical optics 
field over the reference sphere. The results are indicated ±n figures 
17 and 18 as small departures from the heavy dark line^ v;hich represents 
the values used to generate the fouricr coefficients. 

These fourier coefficients vjere then used vzith the series solution ^ 
given by equations 3,16 and 3«17, to evaluate the fields at the reference 
sphere. The trends and accuracy of this series solutioti^ vdien compared 
to the geometrical optics field initially calculated and subsequently 
used to generate the necessary coef ficients , appears to be reasonable 
V7hen viewed in figures 19, It is to be recalled that the series V7as 
truncated at a value of terms much less than infinite in order to avoid 
generation of highly reactive fields. The presentation of figures 19 
involved 390 terms. As the series* accuracy improves \7Xth the number of 
terms, better agreement between the geometrical optics and resulting 
synthesized fields would theoretically be f orthcoraing with continued 
summation. It is to be noticed that the agreement of phase appears 
quite good for the tangential fields and satisfactory for the radial 
component. Similar agreement occurs in amplitude. The basic boundary 
condition that V7as involved required matching of the tangentical fields 
so as to force the field to zero at the reflector surface. Hence, the 
radial component is also determined from this boundary condition, and 
should contain some variation for less than an infinite number of terms. 
The phase of the radial component is similarly determined from the 
real and imaginary parts of the component, and for the radial component 






of near-zero magnitude, the error involved in the difference of small 
numbers begins to override the solution. Hence, for values of radial 
component near zero, rapid phase variations v;ere encountered. Those 
were not plotted as not being significant they added no contribution 
to the v7ork. 
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CHAPTER IV 
CONCLUSIONS 

4.1 As noted by the graphical representation of solutio-n offered in 

figures 19, the technique explored bif this thesis has yielded results 
that indicate its accuracy and applicability, A geometrical optics expan- 
sion from the surface of a spherical reflector to a hypothetical ref~ 
erence sphere can be utilizedto generate necessary coefficients required 
by the series solution of spherical harmonics that constitute the 
expression for the electromagnetic fields in the region of concern. 

With regards to the efficiency of the technique, an investigation of 
various eccentric origins from 7 . - ,5R to ,6R indicated that the radii 
of the ensuing reference spheres v:ere on the order of R' “ .6R for 
the reflector chosen for investigation by this thesis. Thus, evaluation 
of the series solution at the reference sphere for purposes of compar-" 
ison to the initially determined geometrical optics field involved a 
reduction in the number of terms that vrould have been required \7ithout 
the eccentric coordinate system by some 40 %. The reduction in terms is 
due to the direct relationship bctv7een the radius of observation and the 
number of terms required. Due to the flexible nature of the technique, 
choice of eccentric origin to include regions of energy concentration 
near the feed sphere result in the feed radius on the order of 8 X 
for the reflector of 93 X considered by this study. 
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Previous investigation ^ coustrairAcd by a fixed origin and con-- 
centric spheres ^ as given by figure found it necessary to keep the 

radius of the feed sphere on the order of half the rridius of the ref lac 
tor* Thus V7G consider that this technique has provided a substantial 
reduction in the number of terms involved from 



M kr « ( 2 ',t/X )( 93A /2) 



to 



H - kr ( 27 t/X ) 8 X 



indicating a savings of a factor of approximately an order of magnitude 
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4.2 A next logical s,te,p in pursuance of this technique 
would be to concentrate on the series solution at the 
radii of various feed spheres. This study suggests that 
if the fields given by figure 16 can be duplicated over 
the surface of the feed sphere by some source within the 
sphere the desired illumination of uniform plane wave 
over the aperture of the reflector will be realized. The 
calculation of the fields at the aperture by an alternative 
measurement technique, or by experimental methods, if the 
fields over the feed sphere can be realized by some practi- 
cal feed system, would be a most worthwhile endeavor. 

Further study for a series of parameter variations for 
one reflector could, as well, yield a meaningful extension 
of this technique. Based upon the field distribution over 
the feed sphere, design and experimentation with a feed 
system could provide the basis for additional study of this 
topic at various academic levels of concern. 
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/iPPERDIX I 



SIMPLIFICATION AND MAIUPULATIVE DETAILS OF CERTAIN EXPRESSIONS 



D e r :i-v at j on oi: Cq ur'.t.ton 2^12 



Frcn equation 2.11: 



N/p- "Sxvx^ 



A.l 



Taking the dot product with N to both sides of A.l 



l/p = "N «S X V X S 



applying the vector identity 



A*(B X C) -- C.(A X B) 



then 



1/p « - (N X S).(V X S) 



A. 2 



By equation 2.10 



S - VL/n ' 
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Substituting, am! oporating 



V X S V X = V X VL -I- _J_ V X VL 

Ti n n 

But V X VL = 0 GO 



V X S = V X VL • A. 3 

T1 



Substituting A. 3 into A# 2: 



l/p “ “(N X S)*( V X VL) 

T1 



A* A 



We substitute 



V 



1 



n 



n 



l/p = -(N X S).( - J__ V(ln ,1 ) X VL) A. 5 

n 



But by the vector identity 



A X B • C X D (A*C)(B*D) (A*D)(B«C) 
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I 

A 

1 







♦ 



1/p - 



1 



( (M'Vln m)(S‘VL) - (N»VL)(S-Vln n) ) 



A. 6 



n 



Frcr.i 2,105 vjhere S*VL = | VL [ ^ « n 

n 



1 /p “ ( K«Vln n - (N*VL)(S*Vln n) A . 7 

n 



The second terra in A , 7 is identically zero since by A ,1 



N = - p(S >: V X S) 



and from 2.10 



VL «= Sn 



then substituting the above into the second terra; 



(N»VL) =-■ -S*(S X V X S) 
Pn 



Again using the vector identity 



A* (B X C) C‘ (A X B) 



lOZ 




0 



I 



wa GGG 



V X S 0(3 X S) - 0 



TheraforCj from A«7 



1/p - M®Vln n 



Dorivatzion of e quati on 2t37 * 



From equation 2.34 



X - ( z - 2 ) tan20 

o 

( z^- “ 2zz -F* z^ ) tan^*20 
o o 



From equation 2.35 



= x2 ‘h « 2zcR •}- c2r2 



Substitution of A. 9 into A. 10 



R’2 = ( ^2 _ 2zz ‘I* z2 )taii2 20 -J- 
o o 



A. 11 



= (tau220 -f- l)z2 - 2z(z^tan220 - 



A. 8 - 2.12 



A. 9 



A. 10 



“ 2zcR 'I' c2r2 



h cR)i- 

o 
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Mu3-tlplying A* 11 by cog^20 



R'^cos^20 - (sin^?0 -1- cos^20)2^ - 2z (r.^^sin^20 -r cRcog^*20) 



-r 2 :^- Gin^20 -f c^'R^co3^2G 
o 



A.12 



Solving A, 12 for z : 



0 r= 2 ^* - 2z(z sin^20 cRcos^20) -h z ^sin^20 “ cos^20(R’^ « c^R^) 
o o 



From 2,36: 



R 



sin20 



z 

o 



sin0 



z 

o 



sin0 



R 



Gin20 



A.IA 



Substituting A, 14 into A, 13 J 

- 2z (sin0siri2OR + cRcog^20) sin^GR^ - cos^20( R.*^ « c^R^) - 0 

A, 15 



10 ^ 



Ael3 



NorTiia.lizlng A. 15 to the rndiuc of reflectoi: R; 




f 



R’ 

n 



R 



R’ 

R 



- 2z (sin0sixi2G *1- c^coc^2G) -I sin^6 -cos^20 (R.’^ 
n n n 



c^) 0,Ael6^^'2c37 






Normalising equation 2t41, 2.42: 



\7here 



R 



In 



/r 



tm 



2 

cn 



T X, 



mu 



X ) 
cn 



s-2^' 



A. 17 



z ~ 


cog6 




mn 


X = 


slnO 




inn 


z - 


1 ( 


3cos0 - cos3G) 


cn 


A 


X « 


' ( 


3sin0 «• cinSG) 


cn 


A 



z “ sinG 
sin20 



Substitution and expanding /i.l7: 



r2 

In 



( sin 6 
A 



A. 18 



A - 




1 ( 2cos9cos30 -t- 2sln6sin30 + 2 ) 

16 



1 ( 2 + 2cos20 ) 



A. 19 



but cos20 =■ 



then 



16 



1 - 2sin^0 



A ( 1 - siu^O) 
16 



CO5^0 

A 



COG0 A,20~2.A1 



2 



Similarly, by expanding 2.A2 



( COS0 ~ 



sin0 



+ siu^0 



A. 21 



in20 



cos^O -I- sin^O -'r ^inO_ 
siu^20 



2cos 0 sin0 

ciu20 



p2 


= 1 - 1 + sirO 

si.n 20 


Then 







= 1 A.22-* 2.44 

2 CPS0 
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u 
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/vPPEMDIX II 



COMPUTATIONAL TECHNIQUES RELEVMT TO PROBLEM SOLUTION 



As noted in section 2*5 ^ expressions for fields \7cr(i prograrnxaed in 
FORTRAN to yield a nunerical solution to equation 2«52e A cop}’' cf the 
progran utilised follov/s this discussion* Briefly, the lot;ical progression 
is outlined; 

1* For a given reflector of radius R and subtended angle 
On, equation 2*31 is used to calculate the reference 
sphere radius R’ for a particular choice of parancter 
c* 

2. Equation 2* AO is programmed to yield the angle 0’ 

neasured fren the nev; origin at to points on the 

reference sphere corresponding to rays incident upon 
the reflector surface at a given angle 0. A one-to-one 
correspondence is established between the angles. 

3. Equations 2eAl and 2.A2 arc prograrnned along with 
2.33 to yield values for the field anplitu Ic at var- 
ious angles 6’, neasured fron the nev; origin, on the 
surface of the reference sphere. The arplitude expression 
is then given by the inverse distance relationship of 

the forn of equation 2.17. ^ 

loQ 



A. The polarir,atiori fcaturec ccvcreJ by ccraationa 2c 46-2 ,51 
are progriuried to indicate the effect of redefiaiiif, the 
field conponentG ir; the new coerJiaate syotcr,.io 

5* The pliace cxpreooioii ic calculi-ted for tb.e reflector 
under conaidcration by relation 2.30 and 2.33^ and the 
rcnulting phase is chosen relative to the pliase at 
0-0* - 0 degrees. 

6. The iafomation generated by the previous steps is 
collected and output ^ if desired ^ is nade available, 
representing a nuiierical solution to equation 2.52. 

The graphs included in section 2.5 were prepared fron results generated 
as described above. 

It is noted, however, that the actual field on the reference sphere 
is only of passing interest at this point, and that expansion of this 
field into a cosine fourier series, as given by equation 3.14, is of 
paranount interest; 



E reference = 
sphere 
tangential 



CO 

^ 0 - 

n-0 



•f- d E ir'.<{> 
n 



i^)coBJi0 ' 



to9 



A. 23 



The subroutine FORC vas dcvelopad t:c> ger.oratc coef f icicn-ts for the 

fourier series fron the generated field ccnponcats* A nunerical tecluiique 

for evaluating the coefficients for such a series appropriate for 

44 

a periodic set cf discrete values of a function \;as available and vras 
used as the basis for this subroutine. As a chcckj, the fourier series 
v/as sv:i:nedj indicating accuracy to vjitliin a fe\7 percent V7lien viewed 
overall and cerrpared to the points used to generate the coef f icients « 

The tecliViique follows fron consideration of expansion of a periodic 
function f(0) in a feurier series j 



f(e) 



CO 

a /2 -I- Z a cosnG •}* b siiraG 
o , n n 

n=l 



Ac 24 



V7here fron orthogonality it nay easily be shewn that the resviltant 
coefficients arc given by: 



a 

n 



f f(G) cosnO clG 

0 



n^O p 1 j e • c t A. 25 



As the case of interest involves synnetry about 0-0, it is obvious 

that a cosine series \;ill suffice to fully express the function and 

b will be identically zero for all u. 
u 



no 




\ 



The integral exprennion A* 25 can be related to a sun :atiou of diGcrcte 
values by approxi-nations of basic calculus j. v;aere the interval fren 
0 < 0 < 2 tt is divided into 2N intervals ^ providing 2N d- 1 data points c 
Thus V7C can fornulatc a relationship betv/een the nuriber of intervals 
and the spacing involved where 



2N == 360°/ DEL 



DEL bcin^ the spacing bet\;een successive points. Thus we approxinate, 
for a sufficiently large nunber of points; 



do K A0 “ 



2tt 



-2N + 1 



0 = nA0 



2Trn 



v;here n= 0(1.2, ,.,,2il 



2N + 1 



Substitution of the above approxinatious into A, 25 results in 



a 

n 



2 

2N + 1 



2H 

I 

n^O 



f( 



2;t n 



2U + 1 



cosn 



2^n 



2N -i- 1 



A. 26 



IK 



Thus ve Gee the function nunt be nac!e availab.le in Gteps of angle j, 
detcrninecl by the index Uj or at the follov/ing discrete points: 



0 

n 



2iin 

ZN -h 1 



n 0,1,2,*. .2N 



Therefore, the following steps arc outlined to provide an iv\dication 
of the progression follc\7cd through the subroutine: 

1« Since the series requires even spacing betv;ecn the data 
points, it vjas necessary to interpolate between the points 
generated on the reference sphere by successive steps of 
the angle along the reflector surface. As equation 2. AO 
allcv;ed no sinple inverse relationship between the t\;o 
angles 0 and 0’, snail steps of the angle on the reflector 
surface were taken ( spacing - ,2 degrees ) to insure the 
. . accuracy of linear 5.nterpolaticn between the discrete values 

of angle 0* generated. 

2. The resulting field expressions were arranged in coripcnent 



forn. and subsequently expanded fron the exponential represent** 
ation to the trigonometric forn : 



\;Iiere - A( cos\0 *1- jsin\^i ) ~ AH -!- jAI 

3« The va3,ues of field conpoucntc^ tvo real and two 
iuaginary, were then used ,t.o geaerate a conplex 
fouricr coefficient of the forv.i 
C = CR jCI 

by the technique previously described. 

4. The answers of interest were these coefficients which, 
when properly inserted into an existing progran for 
the series solution to field expressions 3.16 and 3.17. 
The fourier ccefficionts c and d are noted in equation 
3.16. 
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